TECHNICAL REPORT R-8 

MASS TRANSFER COOLING NEAR THE STAGNATION POINT 

By LEONARD ROBERTS 

Langley Research Center 
Langley Field, Ya. 


i 




TECHNICAL REPORT R-8 


MASS TRANSFER COOLING NEAR THE STAGNATION POINT 1 


By Leonard Roberts 


SUMMARY 

A simplified analysis is made of mass transfer 
cooling — that is } injection of a foreign gas — near the 
stagnation point for two-dimensional and axisym - 
metric bodies. The reduction in heat transfer is 
given in terms of the properties of the coolant gas 
and it is shown that the heat transfer may be reduced 
considerably by the introduction of a gas having 
appropriate thermal and diffusive properties. The 
mechanism by which heat transfer is reduced is 
discussed. 

INTRODUCTION 

The reduction of heat transfer near the stagna- 
tion point of a blunt body (two-dimensional or 
axisymmetric) is of primary importance when the 
body attains high velocity. One method of cooling 
is to introduce gas of high specific heat into the 
laminar boundary layer near the stagnation point 
and to allow this gas to flow over the nose of the 
body so that large amounts of heat arc convected 
away from the nose. 

A number of exact solutions based on various 
assumptions of fluid properties are available for 
stagnation-point heat- transfer with no injection. 
(A comprehensive list of these references may be 
found in refs. 1, 2, and 3.) It has also been shown 
that the variation of the product of density and 
coefficient of viscosity across the boundary layer 
has an important effect on the heat transfer (ref. 4). 

The boundary-layer equations for a binary 
mixture are well established (ref. 5) and several 
methods are available for the evaluation of the 
thermal and diffusive properties of a binary 
mixture (ref. 6). 


1 Supersedes NACA Technical Note 4391 by Leonard Roberts, 1958. 


Exact solutions (refs. 1 and 3) have recently 
been obtained for air-to-air injection near the 
stagnation point; these solutions have required 
lengthy numerical integration techniques and are 
necessarily limited to discrete values of the param- 
eters involved (for example, the rate of injection 
of mass and the viscosity-law assumption). 

The purpose of this paper is to present a simpli- 
fied anal} r sis by which the effect on the heat 
transfer of injection of air, or a foreign gas, may 
be found and, also, to show more clearly the im- 
portant thermal and diffusive properties of the 
foreign gas if it is to be effective as a coolant. 
Application of the present method to ablation cool- 
ing is made in reference 7. 

SYMBOLS 

x coordinate along wall 

y coordinate normal to wall 

z transformed ^-coordinate 

7 arbitrary value of y outside boundary 

layer 

Z arbitrary value of z outside boundary 

layer 

u component of velocity in ^-direction 

v component of velocity in y-direction 

pv 

V modified velocity component, — 

Pw 

U free-stream velocity in r-direction 

C constant in velocity distribution 

T_ temperature 

T mean temperature (eq. (43)) 

W concentration of foreign gas 


l 



2 


TECHNICAL REPORT R-8 — 'NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


W mean concentration (eq. (42)) 

p density of mixture 

p, coefficient of viscosity 

k thermal conductivity of mixture 

D n coefficient of bina^ diffusivity 

m rate of mass injection per unit area of wall 

c v specific heat at constant pressure 

7*, mean specific heat (eq. (47)) 

N Nu Nusselt number 

Ti Reynolds number 

Np r Prandtl number 

N Sc Schmidt number 

q heat-transfer rate per unit area 

5„ vclocity-boundary-layer thickness 

5 t thermal -boundary-layer thickness 

8 W conccntrat ion-boundary-layer thickness 

K function of Prandtl number and Schmidt 

number (defined by eq. (40)) 

II eff effective heat capacity of foreign gas 

Subscripts: 

e external flow, near stagnation point 

w wall 

c coolant, far inside wall 

0 no injection 

1 foreign gas 

2 air 

b body 

DESCRIPTION OF THE BOUNDARY LAYER 

The flow considered is that shown in figure 1. 
In the steady state the laminar boundary layer 



Figure 1. — Flow configuration. 


near the stagnation point is a thin layer of fluid in 
which the velocity, temperature, and concentra- 
tion of the foreign gas vary rapidly from the 
external stream values to the wall values. In the 
neighborhood of the stagnation point the three 
superimposed boundary layers (that is, the 
velocity, temperature, and concentration bound- 
ary layers) have constant, although different, 
thicknesses. 

It is assumed that the component of velocity 
parallel to the wall is linear in x and that the nor- 
mal component of velocity, the temperature, the 
concentration of the foreign gas, and the prop- 
erties of the mixture are all functions only of y , 
the distance normal to the wall. (This assump- 
tion results directly from the “similarity” nature 
of the flow and the absence of thermal and concen- 
tration gradients along the wall.) It is also 
assumed, consistent with the foregoing variation of 
velocity components, that the coolant gas is 
injected normally at the wall with a velocity 
independent of x. 

When a given amount of coolant gas is injected, 
it diffuses through the boundary layer and is con- 
verted with the air, as a mixture, under the action 
of the pressure gradient imposed by the external 
flow and the shearing stress due to the presence 
of the wall. The concentration of foreign gas at 
the wall is uniquely determined by the rate of 
injection. The extent of the boundary-layer 
shielding depends upon the specific heat of the 
foreign gas and upon the wall temperature as well 
as the coolant and stagnation temperatures. The 
shielding also depends upon the manner in which 
the coolant gas diffuses through the boundary 
layer. The analysis shows the relative importance 
of these quantities. 

ANALYSIS 

The method to be used is as follows: the wall 
conditions for coolant injection are formulated in 
a simple manner (and the results are justified by a 
more detailed consideration of the flow within the 
porous wall presented in the appendix); use is 
then made of previous “exact” results for stag- 
nation-point heat transfer with no injection in 
order to determine the gross characteristics of the 
thermal and viscous boundary layers near the 
stagnation point. These results are then used 
in a simple approximate integral method to 
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determine the effect of injection on the heat trans- 
fer to the wall. 

THE COOLANT FLOW 

The boundary conditions at a porous wall 
through which a coolant is injected may be ob- 
tained very simply by ignoring the presence of the 
solid part of the wall; the justification for this is 
given in the appendix where coolant gas flow within 
the porous wall is considered in more detail. 

When the steady flow of coolant toward the sur- 
face of the wall is considered, the volume taken up 
by the solid wall being neglected, it is seen that 
there is a balance of diffusion and convection 
which governs the flow of mass and heat within 
the wall. 

The transfer of mass is given simply by 


It is important to note here that, even though the 
specific heat of the mixture c p is given by 

c P =Cp' l W+c P ' 2 (l— W) (4) 

(with c Pt i and c Pt2 constant), the value used cor- 
rectly in equation (3) is c Ptl ; in unit time the heat 
transferred at any point y in a direction away from 
the surface is sufficient to raise an amount of 
foreign gas of specific heat c pA through the tem- 
perature range T— T c . 

Evaluating equations (2) and (3) at the wall 
surface gives the following boundary conditions 
which are required for the solution of the bound- 
ary-layer equations: 


pv=p w Vu,=m ( 1 ) 

where rh is constant. 

The diffusion of air inward from the surface is 
balanced by the convection toward the surface: 

~ pDl2 ly = (1 “ W ^ ( 2 ) 

Diffusion of air Convection of air 

from surface toward surface 

since w is the concentration of the foreign gas and 
(1— IT) is that of air. 

Similarly, the transfer of heat is given by 

k~ = c p,i(T—T c )m (3) 

Diffusion of heat Convection of heat 

from surface toward surface 


L)'» (6) 

Detailed study of the flow within the porous wall 
yields the same boundary conditions. (See ap- 
pendix.) 

THE BOUNDARY-LAYER INTEGRAL EQUATIONS 

The integral equations which describe the 
transfer of mass and heat in the boundary layer are 
derived in a simple way without reference to the 
general differential equations. 

Consider a small rectangle of height Y (7" 
greater than any boundary -layer thickness) and 
length Ax } in which p, TC, and T are independent 
of x , near the stagnation point of a two-dimen- 
sional body. (See fig, 2.) The continuity of mass 
may be expressed as follows: 


— p(Y)v(Y)Ax + mAx — 

Flow in from Injection 
external stream 


\{ Y pudy~\ 

LJo Jx+Az 

Flow out at 
x~\~ Ax 


[f>*l 


Flow in 
at x 


This equation, in the limit as Ax->0, becomes 

-p(Y)v(Y) +rh=£ J" P u dy ( 7 ) 


and 


n u 

u—C jj x 


UL 

When the form U=Cx of the external velocity is whcre JJ is independent of ar, then equation (7) 
used and it is noted that reduces to 


- P (Y)v(Y)+m=cJ o r pgdy 




( 8 ) 



4 


TECHNICAL REPORT R8 — “NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Since the density p is a function of concentration 
and temperature, it is convenient to introduce the 
quantities 



dz= 



J 


(9) 


so that equation (8) becomes finally 


— Pu>V(Z)+m— 



( 10 ) 


Similarly, the equation for the transfer of foreign 
gas becomes 


m 

Gas injected 


y 



Figure 2. — Boundary-layer mass balance. 


Cp^Wjjdz ( 11 ) 

Gas convected 


The equation of heat transfer is written 


— C v , 2 ( T e — T w ) P U V{Z) 

Convection of heat 
from external stream 



Transfer of 
heat to wall 


C Pw j\ p (T-T a ) jj dz 

Convection of 

heat in boundary layer 


( 12 ) 


The quantity V(Z ) may be eliminated from equa- 
tion (12) by use of equations (4), (10), and (11); 
thus, 

Cp.f g *C,(T,-T) T u .) rh 

(13) 

Alternatively, by use of the boundary condition 
given by equation (6), equation (13) lias the form 

Cp w f*c,(T,-T) JJ dz=c pA {T'-T c )m (14) 

The equations for transfer of heat and mass 
near the axisyrnmetric stagnation point can be 
derived in a similar manner and arc identical with 
equations (10) to (14) where x is again the dis- 
tance from the stagnation point along the surface 

but 2=2 j - - dy. 
d 0 Pw 

HEAT-TRANSFER RELATIONS FOR NO INJECTION 

The rate of heat transfer is first given in terms 
of Nusselt number N NUtW , Reynolds number /? u ,, 
and Prandtl number N PrtW for which the usual 


definitions near the stagnation point arc used: 


and 


at — 

1 ’ A" u . w ~~ “ 


X 

T 

1 £ W 



Rw= 


Pvd ^ f P w 

Pw Pw 


X 


2 


\r ^ r 

h u 


P.2 


The rate of heat transfer to the wall with no 
injection q Q is now 



(15) 


Exact evaluation of the dimensionless parameter 

('TTI/r) has been made by many investigators 
\ a» /o 

(for example, refs. 1, 2, and 4). The following 
formulas were found to be in good agreement with 
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previous exact results (ref. 1) under the assumption 
of constant pi u: 

For the axisymmetric body, 

(■-%)] <i6) 

and for tlie two-dimensional body, 

(x : ;'). = [ o ' 57O ~ OO65 ('-0]- v '''*"‘ (17) 

Alternative expressions, based on the solutions 
presented in reference 4 in which p\x is not constant, 
are as follows: 

For the axisymmetric body, 

(^7fW) = 0 - 765 (— T (18) 

\ ' So \PuPw/ 

and for the two-dimensional body, 

(t^).-°' 570 (S;)" V ''-- (19) 


by the superscript (1), then equations (16) to (19) 
show that 


/ ATjvh, «>\ / A^vti, mV 1 * 


0.4 


(23) 


This equation shows the effect of the Prandtl 
number on the heat-transfer coefficient. When 
A r Pr , fiJ =l, the thermal and viscous diffusive effects 


are similar so that -^—1; thus, from equation 

0 

j -"Cfr=°(vb)"’ <M) 


The nondimensional thickness <5 Uf0 


(£7 


is now 


assumed to be independent of the Prandtl number 
(since the momentum equation is only weakly de- 
pendent on the Prandtl number) ; thus, when equa- 
tions (23) and (24) are inserted into equation (22), 
the following simple result is obtained: 


BOUNDARY-LAYER THICKNESS WITHOUT INJECTION 

Use is now made of these results (eqs. (16) to 
(19)) to determine the thickness 5 Ui0 of the velocity 
boundary layer without injection. Equation (13) 
with m~0 is written in the form 


1 C z T e ~Tu 
du.oJo T € - T w U 


J_ / MtrV 72 1 /N NUtW \ 

du,0\Pw^J Np Tt to \ H w l/2 /o 


(20) 


and linear profiles arc assumed for velocity and 
temperature in order to evaluate the integral in 
equation (20) ; that is, 


u z 

U 5 Ui0 


T.-T 

T e —T w 


1 





( 21 ) 


where is the thickness of the thermal boundary 
layer. The left-hand side of equation (20) is then 

approximately ~ and equation (20) reduces 

h \o Ut 0/ 

to 

I J / Y /2 * / 

6 \5 Ut0 ) Ko \PwCj Npt.wK K» /2 Jo ' 

If the value of ( w\ieji'N Pr , w =l is signified 


(25) 

d U,G 

This result is analogous to that for a flat plate 
(ref. 8) : 


INJECTION OF AIR 

Consider now the effect of injection of air into 
the boundary layer. The boundary layer be- 
comes thicker because of the increase in mass flow 
which results in reduced gradients in the boundary 
layer. In particular, the heat transfer to the wall 
is reduced since the heat convected parallel to the 
wall in the boundary layer is increased. 

The effect of air injection on the boundary -layer 
thickness is given by considering the mass flow in 
unit time per unit area, as follows: 


Cpw SJu dz 

= 


m 

Mass flow with 


Mass flow with 

Mass 

injection 


no injection 

injected 

which gives 

2 S “ 

- 1 A 4- ™ 

(26) 

when linear profiles 

are assumed for ^ 

(that is, 
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U 

U 


— (no injection) and (injection) )• 

8 U ,0 0 M / 


It is assumed that the relation given by equa- 
tion (25) is still true since the ratio of thermal- 
boundary-layer thickness to velocity-boundary- 
layer thickness depends on diffusive processes, 
that is, on the Prandtl number. 

With c Ptl =c Pt 2 equation (13) may be written 
by use of equation (15) in the form 


1 C 7 T-T u 

5 Jo T^T w U aZ 


v /2 r * w 

&u\PwC) L N Pr ,„RJt*' 


m i 

{p^ w cW\ 


Linear profiles for temperature and velocity 
(eqs. (21)) and the expressions for 8 U (eq. (26)) 
and 5 W , 0 (eq. (24)) arc used to simplify this 
equation, with the following result: 


1 N*u ,1 r 1 /^Nu,ip\ 

n^7 w rj I2 ~n^A rj 12 )* 

1 ^ AT -0.6^ ^ 

V 1 3 Apr '* Jip^C) 112 


(27) 


In terms of the rate of heat transfer q , equation 
(27) has the form 

q = 2 o - (l-lNr r ,.- 0A )c P AT-T K )m 

m 

Heat trans- Boundary-layer shielding by 

fer for no convection 

injection 

The quantity (j Jj Av r)tt r 0,6 ^ (T e — T w ) may be 

interpreted as the average temperature rise in the 
boundary layer of the mass introduced at the wall 
surface. An alternative expression of this heat 
balance, obtained by using the boundary condition 
given by equation (6), is 

qo=c„, 2 [(TV- 7V)+(l N Pr , 6 )(r T w )]m (29) 


When q 0 is expressed in terms of dimensionless 
quantities by use of equations (15) and (23), 
equation (29) becomes 


T-T w 

T~T C 


/V °- 6 m 

_ rr - w (p^«cyv 

/ w\ <V I 1 _ ^ 

V IW 2 A .3 (p*p w O V2 


(30) 


Equation (30) may be used to determine the 


rate of injection required to maintain the wall 
temperature T w at a given desired value when 
T e and T c are specified. 

Before proceeding to considerations of injection 
of foreign gas, the results given by equation (27) 
arc compared with available exact solutions as a 
check on the validity of the several simplifying 
assumptions that have been made. 

Firstly, under the assumption of constant pp, 
equations (16) (eq. (17) for the two-dimensional 


flow) and (27) are used to give as a function 

of the dimensionless rate of mass injection 
- ? “TnTTT A comparison of these results with the 

\PwPw'~ / ) 

results of references 1 and 3 shows very good 
agreement except for the extreme rates of injec- 
tion with Prandtl number equal to unity. (See 
figs. 3, 4, and 5.) 

Secondly, for variable pp, equations (18) (or 
eq. (19)) and (27) are used and the results arc 
compared with those of references 3 and 4. For 
a Prandtl number of 0.71 and no injection reference 
4 gives, for the axisymmetric case, 




4 


and the present report gives (from eq. (18)) 


A A r u ,u 

7?„ 1/2 


•=0.667 



A comparison of the results obtained b} r use of 
equation (18) with the results of reference 3 in 
which the Sutherland viscosity law was used is 
presented in figure 6. 

The effects of injection are also compared when 
PH is variable; again good agreement is found 
except for the high injection rates. (See fig. 5.) 

With this justification of the assumptions that 
have been made, the method is now extended to 
take account of the additional effects which 
result from the injection of a foreign gas (that is, 
gas having properties different from those of air). 


INJECTION OF FOREIGN GAS 

The mechanism by which air diffuses through 
the boundary layer in a binary mixture due to the 
concentration gradient is exactly analogous to 
that by which heat diffuses in the absence of 
viscous dissipation for air-to-air injection. In 
this analogy the coefficient of binary diffusivity 
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A 


\ 


\ 


I 


.5 


ct 


\ 


\ 


Q> 

E 

D 

Y_ 

D 

CL 

\_ 

Q> 

<+- 

</) 

C 

o 


D 

Q> 

X 


.3 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


\ 


p Ref. I v 

Present report (eq. (16)) p 

i 01 , \ 


0 .5 1.0 

Rate-of- injection parameier, 


1.5 


m 


[p w \x w C) 


1/2 


(a) Axisymtnetric case. Np r .u> — 0.7; 


0.05. 


Figure 4. — Effect of air-to-air injection on heat-transfer 
parameter. 


and compared with the heat-transfer equation for 
air-to-air injection, 


f z T-T u m 

P "J,, T-T W U T-T v 

T-T c 


(32) 


The boundary condition at the wall is 


The analogous quantities are presented in the 
following table: 


Binary 

diffusion 

Thermal 

diffusion 

TT 

6- E> 

1 1 

1 TT 

S -3 

1 i 

ir* 

T e — T w 

T t - T c 

/>» 

k 

pCr>,2 

A r S c 

Npr 


It may be verified that all the boundary condi- 
tions at the wall, in the external flow, and at large 
distances into the interior of the wall agree with 
the foregoing analogy. 

Because of the foregoing comparison the fol- 
lowing results may he written immediately: from 
equation (25), 

(35) 

°u 

where 5 W is the thickness of the concentration 
boundary layer, and from equation (30) 


W w = 


AT 0.8 ra 




”+1 


m 


36) 


3 (pu^C)^ 

/V \ (U 

The value of \jfvr) to be used in c q iiallon 


(36) is necessarily that given by equation (18) or 
equation (19) since pii is now a variable quantity 
and depends on the concentration of foreign gas 
in the boundary layer. 

When the gas injected is different from air, 
the amount of shielding by r convection in the 
boundary layer depends critically upon the spe- 
cific heat of the gas c Pi i. 

The heat-transfer equation (13) is rewritten, hy 


use of equations (4) and (6), as 


which is compared with 


[' 


k 


pCp.i 



(34) 



-c p , 2 )W](T-T) ~dz 


=c pA rii{T t —T c ) (37) 
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Rate -of -injection parameter, — Rate -of- injection parameter, — 

{p w fi w C)^^ 


(b) Two-dimensional case. AVr,w=0.7. 
Fioure 4. — Concluded. 


and the linear profile 


WvT 5 W 

is assumed. 

Substitution of equations (21), equation (38), and the definition of N Nu 
results similar in form to those for air-to-air injection (eq. (27)): 


(38) 


into equation (37) gives 


1 Aa 


1 


A 7 


Pt , W 


R ' /2 Nr 




6 ) 


m 




(39) 


where 


K = 6 


] f z W T-Tu J 

a Jo w a T-T.C 


(40) 
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Rate -of -inject ion parameter, — 

(p^C) w2 

Figure 5. — Effect of air-to-air injection on heat-transfer 
parameter for two-dimensional case. Np r . w ' -- 0.7; 



Fig 17 re 6. — Variation of 'Z 1 !',! with - Axi symmetric 

/T u? /J PnVw 

stagnation point; no injection. 


and is a function of N FTiW and N SCtW . 

Equation (39) may be expressed in tlie alternative form 


5= So— 


The second term on the right-hand side of equa- In terms of these moan values equation (41) 

tion (41) has the form becomes 


c P.r 


l-KN Sc , w °* 

a, -0 ' 


l-^A r /v.»-° 




(41) 


c P JV+c P Al-W)](T-T w )n 


m 


where W, the effective mean concentration of the 
foreign gas in the boundary layer, is 


i TSA7 0,f. 

1 IV± \ g Ct w 


(42) 


1— JAW 


and where T, the mean temperature in the bound- 
ary layer, is 

r=r w +(i-|AVr, (4.3) 


q = <7o — [C/mTJ +^^, 2(1 lb )](T 7 T w )rh 

Heat Heat transfer Boundary-layer shielding 

transfer for no 

injection 

(44) 

Thus, in order to achieve the maximum amount 
of shielding for given temperature conditions 
and mass injection, c Pii should be made as largo 
as possible and, in addition, T should be made as 
large as possible; if c Pi C>c l}t 2 then IF should be 
large, and if c Vt \<ic p ^ then W should be small. 
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T—T 

The variations with N Pr w and N Sc ,w of ^ — sr 

< J- e 1 w 

and W are shown in figures 7 and 8, respectively. 
It is seen that W increases as N Se>w decreases and 

Y—T 

that increases as N PftW increases. 

e to 

Maximum shielding in the boundary layer is 
therefore achieved when, for c Vi ^>c p ^ N PTtW is 

i.O r 



o A .8 1.2 1.6 

Prandtl number, N Prw 
f- T 

Figure 7. — Variation of rp with A T p r .w' 



Figure 8. — Variation of W with Nsc.w {W is virtually 
independent of N Pr , w .) 


large and N Sc>w is small and, for c Pr i<c p , 2 , A T Pr ,u> 
is large and A r Sc ,w is large. 

These diffusive effects on the shielding are 
explained qualitatively in the following way (fig. 
9): the convective shielding is most efficient when 
the gas of higher specific heat is transported in 
the regions of highest temperature and velocity, 
that is, in the part of the boundary layer farthest 
from the wall. 

Equation (25) shows that when N PrtW is large 
the velocity boundary layer is thicker than the 
thermal boundary layer and thus more of the 
hot gas mixture is convected. When c Pt i^>c Pt2 
and N Sc ,w is small, equation (35) shows that the 
concentration -boundary-layer thickness is greater 
than the velocity-boundary-layer thickness — that 
is, the foreign gas of higher specific heat diffuses 
quickly through the boundary layer before being 
convected. On the. other hand, when c Ptl <^c Pi 2 
and Nsc.w is large, the foreign gas remains near 
the wall and displaces the air of higher specific 
heat into the hot fast region where it is convected. 

An important consideration in any cooling sys- 
tem is the weight penalty which is incurred in 
order to maintain the wall at a desirable tempera- 
ture. The relation between the rate of mass flow 
m and the temperat ures T C} T W) and T c is obtained 
by substitution of the expressions 


and 


q=c Pfl (T w — T c )m 

q^c„ i2 (T'-r w ) (p„M,X r y /2 ~ 

J-v Pr, tc \ n w so 


(from eqs. (G) and (13), respectively) into equa- 
tion (41). This substitution results in the follow- 
ing expressions ^when either equation (18) or 

equation (19) is used for 
For axisymmetric flow, 

0.765 (j^'T -Vpr, w _0 ‘®(Ptrf‘»C) t/2 c p , 2 (T e — T w ) 

\Prcl*w/ 


m- 


Cp, \{T W — T c )+c p (T — T w ) 


and for two-dimensional flow, 


(45) 


0.570 (^Y' 4 A T P r.„-™{p v p w O)'i*c v , 2 {T t - T m ) 

■_ \PwPw/ 


m= 


c v .i{T a —T e )+Cj,(T—T K ) 


( 46 ) 
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(a) c p , 2 >c J>( ]; maximum shielding when 

(b) c p ,i>Cp, 2 ; maximum shielding when 

Figure 9.-- Qualitative description 



Np r ,w is large and A T sc,w is small, 
of the shielding mechanism. 


where __ _ 

c v =c v ,iW-\-c Vti {\ — IT) (47) 

Probably the most important parameter is the 
total effective heat capacity 77 e// of the coolant; 
this parameter is a measure of the total heat 
absorbed when unit mass of the coolant is used 
and is defined as 

H eff =^=c Pjl (T a - T c ) + c P (T- TJ (48) 


eq. (45) or eq. (46)) on the rale of mass injection; 
both these quantities depend on the concentration 
TTV (See eqs. (47), (42), and (35).) 

The following method, however, gives results 
with relatively little computation: 

(1) The parameters N Pr ,w and N Sc ,w are found 
in terms of W w for the particular binary mixture 
under consideration. (Simple empirical methods 
are given in ref. 6.) 

(2) Equation (43) 


METHOD OF CALCULATION OF RATE OF MASS 
INJECTION 

The problem of most interest from practical 
considerations is that of determining the rate of 
mass injection m necessary to achieve a desired 
wall temperature T w when the stagnation temper- 
ature T e and the coolant temperature T c are given. 
This calculation is made complicated by the de- 
pendence of the mean specific heat c p (required for 


T T 1 

1 1 * ^ 1—1 A 7 - 0.6 

T,—T w 3 >ir ’” 


is used to give T (fig. 7). 
(3) Equation (42) 


W- 


1 -KN SCt ™ 
l-\N Pr , a -™ 
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is used to give W (fig. 8). 

(4) Equation (47) 

c v — c p . 1 + c p, a ( 1 ~ ) 

then gives c p (fig. 10). 

now be made from the results of the substitution 
of equations (43), (42), and (47) into either equa- 
tion (45) or equation (46). 

(6) A second plot of - ) against 

(PwHwCy^X PeHeS 

W w can be obtained as a result of converting 
equation (36) into the following form for the 
axisymmetric case: 

m / 0.76511^ ^ 

(Ptcpw^) 1 ^ 2 \ Pepe / AT 0.6 ^ TIT 

1 T Sc, v> — g ** w 


or into the following form for the two-dimensional 
case: 


m 


{pwVwCy 


/ P wP uT\ 
! \ PePc ) 


0.570W. 


A 7 . 




(50) 


The intersection of the curves found in steps 



Figure 10.— Variation of — with W for various values of 

Cp, 2 


c p , 1 ( 
C p ,2 


(5) and («) give IT. and 

(7) Now, p w p w can be found from W„ and 7^, 
and p e fjL e from external conditions; thus, m can be 
calculated. 

Step (7) may be omitted if the approximation 


f PwPic^ * 

^ PwPic\ 

\ PePc J 

V PePc ) 


is made, in which case ^t 7 2 is given directly by 
steps (5) and (6). 

In order to see the main effect of — on the 

Cp, 2 

mass-flow requirements for given temperature 
T — T 

parameter jr — , values of the Prandtl number 

J- e * w 

and Schmidt number equal to unity are taken. 
Thus, from equation (43) 

T— T w 2 

T e —T w 3 


from equation (40) and from equation (42) 

TF=|. Using this value for W in equation (47) 
gives 

_3,1 

c v i 4^ p ’ 2 

Then, equation (45) for axisymmetric flow 
becomes 


m 

(PePeO^ 


m / PwPtr^l Y' 4 

(PwPwC) 1/2 \ PePe ) 


0.765 


Cpj/Trc-Tc , 1\ , 1 

C'ATe-T** 2/ t 6 


and equation (46) for two-dimensional flow 
becomes 


m 0.570 

(p^o ,/r (M /ZW 41 V 1 

Cp, 2 vn %) o 


The results are shown in figure 11. It is seen 
that the required mass flow is reduced by a large 

factor for the higher values of — . This is 

Cp, 2 

especially true for the higher heating rates 

( T — T \ 

as t jr — TrT — ^0 ) inasmuch as the shielding by heat 

1 6 1 id / 
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Figure 11. — Variation of mass-flow parameter with 

rp rp 

- w — - rp c for various values of 

1 e -f u? Cp,2 


convection in the boundary layer is large. 

CONCLUDING REMARKS 
An approximate analysis has been presented 
whereby the reduction in heat transfer near the 
stagnation point may be calculated with little 
difficulty when the coolant properties are known. 
The agreement with available exact solutions for 
air-to-air injection is extremely good in view of the 
simple approach employed, and the qualitative 
trends of the results for foreign -gas injection in 
explaining the shielding mechanism suggest that 
the approximate analysis will generally give 
reliable results. It is expected that the simplified 



Temperature ratio, — — — 
' e~~'w 

(b) Two-dimensional flow. 


Figure 11. — Concluded. 


calculation which shows the dependence of the 
coolant mass requirement on wall temperature 
and coolant specific heat will provide a good esti- 
mate for engineering purposes. 

The conclusions of the analysis arc summarized 
briefly as follows: 

Maximum boundary-layer shielding is achieved 
when the gas of higher specific heat is convcctcd 
in the hot fast-moving part of the boundary layer 
farthest from the wall; this requires that the 
coolant gas introduced have 

(a) High specific heat compared with that of air 

(b) High binary diffusivity (that is, small 
Schmidt number) 

(c) Large Prandtl number. 

Of these three properties the first is the most 
important, whereas the diffusive effects are of 
secondary importance. 

Langley Research Center, 

National Aeronautics And Space Administration, 
Langley Field, Va. , July 3 , 1058. 


APPENDIX 


FLOW WITHIN THE POROUS WALL 


For the purpose of the present analysis it is 
assumed that the wall is uniformly porous and 
presents, to the stream, an area A of solid material 
and an area B across which coolant gas is trans- 
ported. 

The manner in which gas is transported through 
the porous material is given by 

- p D 12 ~B + pcTT'B = m(A+B) (Al) 

Diffusion Convection Total gas 

toward surface toward transported 

surface 

since m is the rate of mass introduced per unit 
area of body surface. 

The air contained in the mixture is transported 
according to the following equation: 

- pl) u ~ B = p»(l-W)B (A2) 


may be written (by using eq. (A3)) as 

p 7), 2 

ay 


(A6) 


where pD l2 is a mean value of 

Equations (A5) and (A6) have solutions (similar 
in nature) of the form 


T=T c +(T 9 -T c )e 


k 


(A7) 


which satisfies the condition T— >T C as y— >— oo, 
and 


TF= 1 — (1 —W w )e 


(AS) 


A characteristic thickness 6 r may be defined as 


° t ~L 


T-T c 

T w -T c 


Jy= 


mc„ 


(A 9) 


Diffusion of air Convection 

from surface toward surface 

since there is no net transport of air into the body 
in the steady state. From equations (Al) and 
(A2), 

pvB m ( A + B ) (A3) 

The equation for transfer of beat within the 

c p , i {T— T c ) m (A +B) 

(A4) 

Rate heat absorbed 
by coolant 

Equation (A4) may be written 


porous wall is 

(IT 

(hA+JcB) ~ = 

dy 

Conduction into solid 
and coolant 


k d ^=c pA (T-T c )m 


(A5) 


where k is a mean value of the effective thermal 
conductivity - ^ , ^ — Similarly, equation (A2) 


504995 60 - 


in which the temperature difference T—T c changes 
by an amount e~'(T w — T c ). 

Similarly, the concentration W changes from 
W w at the surface to the value 1 in a layer of 
thickness 0 W (within the porous wall) where 





1—W 
1 —W w 


dy 


pDn 

771 


(A10) 


It is seen that both thicknesses depend on the ratio 


B 

A’ 


It is desirable to limit the relatively high 


temperature region within the wall to a narrow 
region near the surface; the specific heat c p must 
be large and k must be small, that is, both k b and 
k should be small. Considerations of the bound- 
ary-layer flow outside the wall also lead to the 
conclusion that the specific heat c P(I should be 
large and k should be small (so that the Prandtl 
number may be as large as possible). Thus, 
these requirements result in two desirable features: 
(a) Confinement of the high temperatures 
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within the body to a narrow region near the sur- 
face and 

(b) Maximum convection of heat away from 
the stagnation point (as shown in the body of the 
report) 

The conditions at the surface as it is approached 
from within the wall are given as follows: 

From equations (A2) and (A3) 

_(pA 2 f) B=(l— TU_)m(A+B) (All) 

and from equation (A4) 

(**A+*B) (J) K _=^ 7V) m(A+B) 

(A 12) 

where the subscript id — denotes that the surface 
is approached through negative values of y. The 
rate of heat transfer from the stream boundary 

layer is (A+B) so that the boundary 

condition (eq. (Al2)) becomes 

( k S)j =Cp ' l(Tu ~ Tc)7h (A13) 

(1 T 

Thus, there is a discontinuity in -r~ at the surface 


when k w 9 ^k l> . 

Similarly, the condition of no net transfer of air 
to the body is given by 

— (^pl)\2 (A+B) — (1— W w )p w v w (A-\- B) 

(A 14) 

obtained when the surface is approached through 
positive values of y. Since 

P W — V to— ■ == ^ P uVtc 

equation (Al4) may be written as 

-(pl) i2 ^-) = (1 -Tig to (A 1 51 

Tt is noted that v is discontinuous since there is a 
discontinuity in the area available for convection 

given by the ratio The concentration 

gradient is also discontinuous for the same 

reason; this is expected by analogy with the tern- 
d T 

perature gradient Equations (A 13) and (A 15) 

are seen to be identical with equations (3) and (2) 
when evaluated at the wall. 
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